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Abstract. - In connection with recent discussion of topological order and topological phase transi- 
tions in quantum systems, we reexamine circumstances that lead to the appearance of a topological 
glass in certain classical lattice spin models. Local bonding enforces constraints on low energy 
states which organize themselves into topologically distinct classes that break ergodicity but not 
any apparent symmetry as in the usual Landau theory of phase transitions. Various properties of 
such a topological glass are demonstrated using two classical Ising-like models. 



Quantum field theories with non-trivial topological structures have attracted a great 
deal of interest in recent years [1-5]. Although the fundamental variables of these systems 
are topologically simple (e.g., ordinary spins or boson/fermion fields on a periodic lattice), 
their low energy properties are characterized by a non-trivial topological order and order 
parameters that are non-local functions of the original microscopic variables. The effective 
variables may take the form of fermions or bosons dressed by Chern-Simon field [1] or string- 
like objects [2,3]. Topological excitations with fractionalized charge (or spin) quantum 
number often appear in these systems as a result of the topological order. The study of this 
class of problems is important for a better understanding of some of the more exotic states 
of matter such as the FQH states. 

Many existing examples of topological order are from quantum systems with a non-trivial 
Berry phase structure. Non-trivial topology may also arise, however, from the spatial bond- 
ing pattern in the ground state which can be understood from a purely classical viewpoint. 
Using simple and clear examples from statistical mechanics, we examine in this paper how 
various aspects of topological order, such as ground state degeneracy, reduced phase space, 
topological excitations and topological phase transitions, come together in such a scenario. 
Whereas topological order is often characterized by topological excitations with fractional- 
ized charge (or spin) quantum number in quantum systems, it manifests itself as breaking 
of ergodicity in classical systems. These related but also distinct features are noteworthy as 
they contribute to a more complete identification and classification of topologically ordered 
many-body quantum states. 

The classical systems we consider have an extensive ground state entropy and no apparent 
spatial order in the usual sense of the word at any temperature. Nevertheless, the low- 
energy part of the phase space decomposes into topologically distinct sectors, leading to the 
breaking of ergodicity at sufficiently low temperatures. In this respect, they resemble the 
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glass state in disordered systems, though the energetics and dynamics of excitations that 
restore ergodicity are expected to be quite different in the two cases. This point will also be 
discussed in some detail below. 

The AKLT state. — It is easiest to illustrate the concept of a topological glass with 
concrete Hamiltonians. We first consider an example which connects topological glass states 
in classical systems and topological order in quantum systems. The example is a classical 
analogue of the Afflcck-Kcnncdy-Licb-Tasaki (AKLT) Valence Bond Solid state [6] in quan- 
tum spin chains and is defined by the following classical one-dimensional Hamiltonian, 

Hs = J ^ CTiCTi+l, (1) 

z— odd integer 

where J > and cTj = ±1/2 are classical Ising spins (see Fig. 1). Since the Hamiltonian 
is a sum of disjoint terms, (ui) = and there is no long-distance correlation in the system 
at any temperature. Now let us introduce the composite spin variables Sj = oij -\- <T2j+i 
which by construction have possible values Sj = —1, 0, 1. There is no long-range correlation 
between the Sj variables either since they are just local sums of the o-j's. However, a hidden 
string order emerges in the system as the temperature T — > 0. This can be seen from the 
string variable 

N 2N+1 

LmN = ^ Si = ^ CTi. (2) 

i=M i=2M 

At T = 0, the nearest neighbor spins a"2j-i and a2j form ±1/2 pairs and sum to zero except 
those on the first site 2M and the last site 2N ± 1. Therefore the only possible values 
Lmn can take are Lmn = —1,0,1, independent of the site indices N and M. To satisiy 
this constraint, the corresponding Sj variables must display the following string topological 
order: starting with an allowed chain of variables {Sj}, we create another (shorter) chain 
by removing all the Sj^s with Sj = 0. The remaining chain, which is an Ising chain with 
Sj = ±1 only, displays perfect antiferromagnetic Ising order [6]. The string order is a direct 
consequence of the non-trivial topology of the phase space for the Sj variables at T = 0. For 
an open chain of 2N ± 2 sites from i = to i = 2N ± 1, the allowed states, in terms of the 
Sj^s, fall into four disconnected topological sectors labelled by the boundary spin variables 
Co, o'2N+\- States in a given sector can be converted to each other through the spin exchange 
within each pair of sites {2j — l,2j), but conversion between states in different topological 
sectors must involve the flip of boundary spins. The string topological order is responsible 
for the appearance of fractionalized s = ±1/2 Ising spins at the two ends of the chain when 
{Sj} are considered to be the fundamental variables. 
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Fig. 1: The disjoint AKLT chain. Horizontal bonds indicate the nearest neighbor coupling of Eq. 
(!)■ 

The antiferromagnetic Ising model on a triangular lattice. — We next consider 

a more realistic and non-trivial example of the topological glass state: the classical anti- 
ferromagnetic Ising model on a two-dimensional (2D) triangular lattice [7]. The model is 
deflned by the Hamiltonian 

H = jY,SiSj + Y,h'{Si;S), (3) 
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where Si = ±1 and denotes nearest neighbor lattice site pairs. The second term 

h'{Si; S) on the right-hand-side denotes all other plausible short range interactions such as 
the second and third nearest neighbor interactions, magnetic field, etc. The model has been 
well studied by many authors [7-13]. We shall first review briefly properties of the ground 
state at h' = 0, and then examine how topology dictates properties of the system at low 
temperatures with J » h', ksT. 

Ground state and topology of the reduced phase space. Due to the geomet- 
rical frustration associated with the triangular lattice, the ground state of (3) at ft' = is 
not unique. Instead, for three spins on a triangle, the energy is minimized if the spins are 

not all parallel and there arc 2x3 = 6 degenerate configurations. Extending the analysis 
to the whole lattice, we see that any spin configuration that satisfies the constraint 

E ^^ = ±1 (4) 

i=l,2,3 

on every triangle in the lattice minimizes the energy. Equation (4) defines a reduced phase 
space for the spin configurations. Since the number of such configurations grows exponen- 
tially with the system size, there is an extensive ground state degeneracy in this model [7] . 

The constraint (4) encodes a topological order which can be illustrated in a number of 
equivalent ways, including the tiling representation, the string representation, the Solid-On- 
Solid (SOS) representation [8] and the dimer representation on the dual lattice [9]. The 
tiling, string and SOS representations are sketched in Fig. 2. The number of rhombuses in 
each of the three orientations, or the number and orientation of the unbroken strings, are 
topological properties of a given spin configuration which cannot be changed via single spin 
flips or spin flips in a localized region. These properties can be equivalently stated in the 
SOS representation described below which is more convenient for analytical treatment. 




(a) (b) 

Fig. 2: Tiling, string and SOS representation, (a) Connecting neigliboring spins that are antiparallel 
results in the tiling representation (solid lines). The string representation is given by connecting the 

midpoints of the horizontal edges of the rhombuses (dashed lines), (b) The directions Sij of nearest 
neighbor bonds in the triangular lattice used to define height increments in the SOS representation. 

In the SOS representation, a "height" variable hi is assigned to every lattice site i 
according to the following rule: first we assign a direction to all nearest neighbor bonds 
such that they point in counter-clockwise (clockwise) directions when travelling through up 
(down) triangles in the lattice [see Fig. 1(b)]. Starting from a site i = with spin Sq 
and (arbitrarily assigned) height Hq, we travel to the next site i = 1 along sqi- A height 
hi = /lo — 2 is assign to site 1 if the spin on the site is equal to So, otherwise hi = ho + 1- 
The construction is then repeated to the third site. Since Eq. (4) enforces two antiparallel 
and one parallel pairs of spins on every triangle, a consistent way of assigning the heights 
is obtained. The procedure can be continued to all sites on the lattice to yield a consistent 
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and unique set of heights from a given spin configuration. Mathematically, the height hi at 
any given site i is related to ho by 



hi = ho-^ '^{ijj-i.Sjj-i) {3SjSj-i + 1) 



(5a) 



where L denotes a chosen path on the lattice from site to z and j = l,2,...,i are sites on 
the path L. Ijj-i is a unit vector pointing from site j — 1 to site j. Notice that the choice of 
path L is completely arbitrary. Conversely, a unique spin configuration can be constructed 
from a given height pattern starting from a given spin 5*0 at site 0. The spin at site i is 
related to Sq by 



(5b) 



Equation (5) defines a one-to-one mapping between {Si} and {/ij} apart from an overall 
height shift /iq- 

To show that the reduced phase space breaks into topologically distinct sectors, let us 
consider the model on a torus (i.e., periodic boundary conditions). For a system of size 
Lx X Ly, where x is along one of the Sij directions [Fig. 2(b)], the boundary condition is 
Sx,Ly = Sxfi and SL^,y = So,y In terms of the height variables, however, the boundary 
condition takes the form hx,Ly = hx,o + 2n and hL^,y = ho,y — 2Lx + 3m, where n and m 
are integers determined by the spin configuration and lie in a certain range set by Ly and 
L^. Spin configurations with diflEerent values of (n, m) belong to different topological sectors 
because any set of local spin flips consistent with (4) do not modify height increments across 
the system. 

Topological excitations and topological phase transitions. — We next consider 
low temperatures ksT <C J. In this case, excitations with all three spins on a triangle 

pointing at the same direction are allowed with an excitation energy 2J. This excitation 
is called a vortex and has a non-trivial topological structure. In the SOS representation, 
presence of a vortex makes it impossible to assign heights /i^'s in a consistent way. A 
discontinuity Ah = 6q is found for the hi variables on the triangle enclosing the vortex 
[q = -|-(— )1 for vortices located on up(down) triangles = charge of vortex]. A line of 
discontinuity (i.e., branch cut), extending from the center of the vortex to infinity or to 
the center of another vortex with opposite charge, is needed to restore a unique height 
configuration (see Fig. 3 for an example). Notice that e'^'^ = 1 and there is no uncertainty 
in the corresponding spin assignment. The topological nature of vortex excitation can also be 
seen in the string representation, which represents the point of annihilation of two strings [12] 
resulting in a difference in incoming string number (from top of lattice) and outgoing string 
number (to bottom of lattice) by two (see Fig. 3). 




Fig. 3: A vortex in tiling, string and SOS representation (an elementary triangle with identical spins 
at three corners). A singular line (dotted line) extends from the center of the vortex to infinity (or 
another vortex). The heights at the two sides of the singular line have a discontinuity Ah = 6. 
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Despite the extensive ground state entropy and the discrete nature of the spins, the 
topological order resulting from the constraint (4) is quite analogous to that of the XY 
model defined by a phase field 6{x) on the lattice. In fact, one may relate the local phase 
gradient V0 in the XY model to the local height gradient for neighboring sites, 

= -^hj-ii^SjSj_i + 1). (6) 

In terms of Ajj-i, the flux through a closed contour L on the lattice is given by, 

Fl=Y1 kj-i ■ ^3,3-1 = 6?, (7) 

where (7 = 0, ±1, ... is the net number of vortices enclosed by L. The reduced phase space 
at T = has Fj, = for any L. A more relaxed definition of topological order which applies 
also at nonzero temperatures is to demand ~ for sufficiently large L, i.e., the vortices 
are in the confined phase. 

Using the XY model as an analogy, we can now understand the low-energy, large-distance 
properties of (3) and the nature of the transition to the disordered phase in terms of vortex 
excitations that form a dilute gas at low temperatures. The statistics of vortices is described 
by an effective Coulomb gas model [11], 

= \Y^ liljOin - fj) +nY^qf, (8) 

i¥=3 i 

where G(r) ~ ,gln(r/a) (a = lattice spacing), qi is the vortex charge, and /x ~ —2J/kBT. 
For h! = 0, it is found that [11, 12] the ground state entropy yields g — 2 which is less than 
the critical value gc = 4 for vortex confinement in the zero fugacity limit. In this case, the 
vortex gas is in the disordered phase at any finite temperature T > (/z finite). A vortex 
confined phase can be stabilized by introducing a weak next-nearest neighbor ferromagnetic 
coupling h' = —JiJ2(ij)eiLnn^i'^3 which reduces spin (or height) fiuctuations and pushes 
g to much higher values as T ^ 0. Consequently, a Kosterlitz-Thouless transition to the 
topologically ordered phase takes place at fesTxT — 8.48Ji <C J, followed by a second 
transition into a (ordered) flat phase at a lower temperature ksTn ~ 4.4:2 Ji [14]. Spin- 
spin correlation changes from an exponential decay to a power-law decay at Tkt without 
breaking of symmetry. 

Additional terms in the Hamiltonian h' do not affect the topological nature of the low- 
temperature phase as long as they are smaller than J and Ji such that vortices remain 
confined. In this case h' can be expressed in the SOS representation using Eq. (5b). Its 
main effect is to create different phases in the system [12, 14, 15] characterized by different 
behaviors of the height variables hi, meaning that /i, has replaced the spin variables to 
become an effective local order-parameter. 

Higher dimensions. — The XY-type topological order exhibited by spin configura- 
tions in the reduced phase space defined by Eq. (4) can be extended to higher dimensions. 
The scheme of Eq. (6) assigns local height gradients according to the nearest-neighbor 
bonding situation, together with a sign convention set by the bond direction vectors Si_j. 
A consistent height assignment for all sites can be made if and only if the zero-flux con- 
dition Fl = is satisfied for every closed contour on the lattice. Quite generally, Fl can 
be decomposed into contributions from "elementary plaquettes" of the lattice, which are 
triangles in the 2D example. Therefore, if Fj, = holds for every elementary plaquette, 
then a topological order of the XY-type is established. 

As an example, let us consider the three-state chiral Potts model on a body-centered- 
cubic (bee) lattice defined by the Hamiltonian [16] 

Hp = -Y,SWi + h<Ti+sJ, (9) 
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where CTj = 0, 1, 2, Sa is one of four nearest neighbor vectors pointing from the center to the 
four non-adjacent corners of the cube, and S{a, 6} = 1 if a = 6 (mod 3) and otherwise. The 
elementary plaqucttcs of the bcc lattice are formed by the four vectors linked hcad-to-tail 
in six alternative sequences. Since each such plaquette contains four sites, at least one of 
the four bonds is frustrated under the energy (9). It has been shown in Ref. [16] that the 
ground state of (9) contains precisely one frustrated bond per elementary plaquette and 
an extensive entropy as in the 2D case. Assigning a height difference of 1 for the satisfied 
bonds and -3 for the frustrated bonds along directions defined by the s^'s, we obtain a 
unique height configuration {hi} from a given Potts spin configuration {cri}, apart from an 
overall height shift. The inverse mapping is given by cTj = hi (mod 3). With the help of the 
height variables, various topological sectors of the ground state can be identified in a similar 
way as in the 2D case. 

Summary. — Summarizing, we propose in this paper a new kind of glassy state in 
classical systems called topological glass which are analogies of topological order in quantum 
systems. The existence of these states are illustrated by two explicit examples. Topological 
order can emerge in classical systems at low temperature as a result of the non-trivial 
topology of the reduced phase space. In our example of the AKLT state, topological order 
in the classical model is directly linked to the corresponding quantum state. The topological 
order in the 2D antifcrromagnctic Ising model on the triangular lattice is of the XY-type 
which has high dimensional counterparts. It would be interesting to explore topological 
order associated with other symmetries such as those discussed in Ref. [17]. Our work is 
only the first step towards bridging the study of topological order in quantum and classical 
systems. 

* * * 

The authors acknowledge support from the Rcscarc;li Grants Council of the HKSAR 
through grant CA05/06.SC04. TKN acknowledges the hospitality of the Kavli Institute for 
Theoretical Physics China at which this work started. 

REFERENCES 

[I] S.C. Zhang, Int. J. Mod. Phys. B6, 25 (1992). 

[2] see for example, X.G. Wen, Quantum Field Theory of Many-Body Systems, Oxford University 
Press 2004). 

[3] F. PoUmann, J.J. Botouras, K. Shtengel and P. Fulde, Phys. Rev. Lett. 97, 170407 (2006). 

[4] A. Kitaev, Ann. Phys. 321, 2 (2006). 

[5] X.-Y. Feng, G.-M. Zhang and T. Xiang, Phys. Rev. Lett. 98, 087204 (2007). 

[6] I. Affleck, T. Kennedy, E.H. Lieb abd H. Tasaki, Phys. Rev. Lett. 59, 799 (1987). 

[7] R.M.F. Houtappel, Physica 16, 425 (1950). 

[8] H.W.J. Blotc and H.J. Hilhorst, J. Phys. A 15, L631 (1982). 

[9] F.Y. Wu, Phys. Rev. 168 539 (1968). 

[10] J.V. Jose, L.P. Kadanoff, S. Kirkpatrick and D. Nelson, Phys. Rev. B 16, 1217 (1977); H.J.F. 
Knops, Phys. Rev. Lett. 39, 766 (1977). 

[II] H.W.J. Blote and M.P. Nightingale, Phys. Rev. B 47, 15046 (1993). 

[12] B. Nienhuis, H.J. Hilhorst and H.W.J. Blote, J. Phys. A17, 3559 (1984). 
[13] V. Elscr, J. Phys. A17, 1509 (1984). 

[14] X.-F. Qian and H.W.J. Blote, Phys. Rev. E 70, 036112 (2004). 

[15] S. Fujiki, K. Shutoh, Y. Abe and S. Katsura, J. Phys. Soc. Japan. 52, 1531 (1983); H. 

Takayama, K. Matsumoto, H. Kawahara, and K. Wada, ibid. 52, 2888 (1983); D.P. Landau, 

Phys. Rev. B 27, 5604 (1983). 
[16] B.M. Forrest and L.-H. Tang, Phys. Rev. Lett. 64, 1405 (1990). 

[17] J. Sethna, in 1991 Lectures in Complex Systems, Eds. L. Nagel and D. Stein, Santa Fe Institute 
Studies in the Sciences of Complexity, Proc. Vol. XV, Addison- Wesley, 1992. 



p-6 



